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Abstract. We calculate analytically the probability density P(t m ) of the time t m 
at which a continuous-time Brownian motion (with and without drift) attains its 
maximum before passing through the origin for the first time. We also compute the 
joint probability density P(M,t m ) of the maximum M and t m . In the driftless case, 
we find that P(t m ) has power-law tails: P(t m ) ~ i m 3 ^ 2 for large t m and P(t m ) ~ U,}^ 
for small t m . In presence of a drift towards the origin, P(t m ) decays exponentially for 
large t m . The results from numerical simulations are in excellent agreement with our 
analytical predictions. 
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Introduction 

In this paper, we derive the probability distribution of a random variable associated with 
a Brownian motion, namely the time at which a Brownian motion attains its maximum 
value before it crosses the origin for the first time. This random variable appears quite 
naturally in different problems such as in queueing theory and in the evolution of stock 
prices in finance. 

Let us first consider, for example, a single-server discrete-time queueing process, 
modelled as a simple random walk [fl [2] via: 

where l n is the length of the queue at time n and £ n 's are independent and identically 
distributed random variables each taking values +1 with probability p (signifying the 
arrival of a new customer), —1 with probability q (indicating the departure of an already 
served customer) or with probability (1 — p — q). In the queueing language, this is 
referred to as the Geo/Geo/1 queue [2]. 
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busy period 

Figure 1. Queue with busy period. 
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Given Z , one calls busy period the period at the end of which the queue becomes 
empty for the first time (see Fig. 1): during such a period, the server always has some 
customers to serve. It is then natural to enquire about the time at which the queue is at 
its longest during the busy period. In the random walk model where the queue length 
l n is the position of the walker at time step n, this amounts to investigating the time at 
which the position of the walker (initially positive) is farthest from the origin before it 
crosses the origin for the first time. 

Another area where the same variable appears quite naturally is in the evolution of 
stock prices in finance. The evolution of a stock price S n with time n is often modelled 
by the exponential of a random walk [31 H]. Starting from its initial value So the price 
evolves with time stochastically. A natural question for an agent holding this stock is: 
what is the suitable time for selling this stock? If the stock price goes below a threshold, 
say R, it is too risky to wait any longer. Thus an agent can wait at most up to the 
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time at which the ratio S n /R crosses the level 1 from above. Within this time, the ratio 
S n /R will achieve its maximum at some intermediate time which is clearly the best time 
to sell the stock. Assuming that the random variable l n = \n(S n /R) performs a random 
walk starting from its initial value l = \b.(Sq/R) > 0, one then wants to calculate the 
probability distribution of the time at which the random walk is farthest from the origin 
till its first-passage time through the origin, i.e., till the time at which S n /R = 1, i.e., 
\n(S n /R) = indicating the first-passage through the origin. 

In this paper, we will consider a further simplified case, namely a continuous- 
time Brownian motion as opposed to the discrete-time random walk in the above two 
problems. For a continuous-time Brownian motion we calculate explicitly, using path- 
integral methods, the probability density of the time t m at which a Brownian motion 
x(t) (starting from x(0) = Xq > 0) is farthest from the origin before it crosses the origin 
for the first time. Algorithmically speaking, for each sample of the Brownian motion 
starting at Xq we stop when it crosses the origin for the first time say at time t = tf 
and locate the time < t m < tf at which the Brownian motion achieves its maximum 
value. Note that both tf (the first-passage time) and t m varies from one sample to 
another. We repeat it many times and then construct a histogram of the £ m 's which 
gives its probability density function P(t m ). Even though the discrete-time problem is 
more relevant, we expect the continuous-time result to provide the right asymptotics for 
the discrete problem. As we will see below, the continuous-time problem, though still 
non-trivial, is easier to handle analytically. 

Note that for a Brownian motion or a Brownian bridge over a fixed time interval 
[0,T], the probability density P(t,T) of the time r at which the process attains 
its maximum is well known [5]. For example, for a zero-drift Brownian motion 
over [0, T] starting at the origin, the probability density P(r, T) = ^g{r/T) where 
g(x) = l/[ffy/x(l — x)] for x G [0, 1] [5]. On the other hand, for a Brownian bridge over 
the fixed interval [0, T] and starting at the origin, the probability density of r is uniform, 
P(t, T) = 1/T for < r < T [5J. In contrast, in our case, the Brownian motion is not 
over a fixed time interval, but rather over a variable time interval [0, tf] where the upper 
edge tf is the first-passage time which itself is a random variable [6] and hence varies 
from sample to sample. 

The statistical properties of the functionals (such as the area, the maximum etc.) 
of a Brownian motion or its variants (such as a bridge, excursion, meander etc.) over 
a fixed time interval have many applications in physics, graph theory, computer science 
and they have been studied extensively (for recent reviews on Brownian functionals 
see |7l IE]). In particular, the area under a Brownian excursion or meander has found 
many recent applications in problems as diverse as fluctuating interfaces 0], graph 
enumeration [10], lengths of internal paths in rooted planar trees [8j [11] or cost functions 
in data storage via the "linear probing with hashing" algorithm [SJ 112] . Similarly, the 
statistical properties of functionals of Brownian motion restricted up to its first-passage 
time (usually referred to as 'first-passage functionals') also have various applications, and 
have appeared recently in many different contexts [TJ [151 IS] > including the computation 
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of the time period of oscillation of an undamped particle in a random potential [T4] and 
the determination of the distribution of the lifetime of a comet in the solar system [8|H~3]. 
The probability density of the area swept by an initially positive Brownian motion till 
its first-passage time was computed exactly in [15J, with an application to queuing 
theory. In this paper, our focus is on the random variable t m which, though not quite 
a functional in the strict sense, is an important random variable associated with such a 
Brownian motion restricted up to its first-passage time. 




Figure 2. Maximum before the first-passage through the origin for the zero-drift case. 

In [15], the authors also computed directly the probability density P(M) of the 
maximum M of a Brownian motion (starting at xq > 0) before its first-passage time 
through the origin, via a "backward" Fokker-Planck method and showed that it has 
a power law behavior P(M) = xq/M 2 where M > Xq. In this paper, we extend this 
work using a path decomposition method that allows us to obtain the joint probability 
density P(M,t m ) of the maximum M and the time t m at which the maximum occurs 
before the first-passage time. By integrating over M, we then get the 'marginal' P(t m ), 
i.e. the probability density of t m . We calculate P(t m ) explicitly both for a driftless and 
drifted Brownian motion. We also compare the results of numerical simulations to our 
analytical predictions and find excellent agreement. 

1. Driftless Case 

We consider a continuous-time Brownian motion evolving via dx/dt = £(£), where £(t) 
is a white noise with = and = 5(t — t'). We start by recalling the 

quick derivation of P(M) given in [X5] : 

Let q(x) be the probability that a Brownian particle starting from x G [0, M\ exits 
the interval for the first time through 0, i.e., the probability that the maximum before 
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the first-passage time is less than or equal to M. Writing 0At(Ax) for the distribution 
function of a Brownian displacement in the time interval At, we have: 

q(x) — / q(x + Ax)<f>At(Ax) dAx. (1) 



Expanding q(x + Ax) for small values of Ax, and using the fact that in the absence of 
drift the mean value of Ax is 0, one finds that q satisfies: 

= 0, q(0) = 1, q(M) = 0, (2) 
whose solution is: 

q(x) = i - ^. (3) 

As mentioned in its definition, it can easily be seen that q(x) also corresponds to the 
probability that the maximum before the first-passage time is less than or equal to M; 
therefore, differentiating Eq. [3] with respect to M gives the probability density of M: 

P(M) = JL. (4) 

To compute the joint probability density P(M, t m ) we proceed as follows. We first 
assume that the maximum occurs at t m and then we split the Brownian path into two 
parts (before/after t m , as shown in Fig. 3) and determine the weights of a path's left- 
hand side and right-hand side separately. Note that due to the Markovian property of 
the Brownain path, once the position of the walker is specified at t m , the weights of 
the left and the right parts become completely independent and the total weight is just 
proportional to the product of the weights of the two separate parts. For the left part, 
we have a process that propagates from x at t = to M at t = t m without crossing 
the level M in [0,t m ] (since M is the maximum) and the level (the origin). For the 
right part, the process propagates from the value M at t = t m to at t = tf where 
tf >t m without crossing the level M and the level in between. We need to be careful, 
however, because, as pointed out in [9], a Brownian walker that crosses a given level 
once crosses it infinitely many times immediately after the first crossing. It is therefore 
impossible to enforce the constraint x(t m ) = M and simultaneously forcing the motion 
to stay below M before or after t m (for a lattice walk, one does not have this problem 
since the lattice constant provides a natural cut-off). Following the method used in [9], 
we introduce a cut-off e by imposing x(t m ) = M — e and consider all paths having a 
maximum less than or equal to M and passing through M — e at t = t m . We compute 
their weight and then let e go to eventually. 

On the right side of t — t m : we have to determine the weight of a path that starts at 
M — e and exits for the first time the interval [0, M] through 0. This is given by Eq. [3j 

l(M -e) = j- (5) 

On the left side of t = t m : we use a path integral treatment with the Feynman- 
Kac formula (as in [9]) giving the weight of a path in terms of the propagator 



x 
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M-e), where # = -I J* + V(x) with V(x) a square well having infinite 
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Figure 3. Splitting probabilities 



barriers at x = and x = M and V(x) = for < x < M (the infinite barriers at x = 0, 
x = M enforce the condition that the path can penetrate neither at x = nor at x = M). 
The normalized eigenfunctions of H labelled by the integer n — 1,2,3,... are ip n (x) = 

^J~^ sin i^jf) with the associated eigenvalues E n = n 2 -K 2 /2M 2 . The eigenfunction 
ip n {x) vanishes at both ends x = and x = M of the box. The propagator can be easily 

M-e) = M*o)MM - e) e~ E ^ and 



evaluated in this eigenbasis, ( x 
one gets: 



-Ht n 



M-e 



2 ^ _ /mrx \ . fnn(M-e)\ 
— > sin — — — sin — e ^m 2 " 

n=l 



M 



In the limit when e —± 0, we get to leading order 



-Ht v 



(6) 



(7) 



n=l 



Taking the product of Eqs. Eland [7J we get the total weight of the path, to leading order 
in small e, 



P(M,t m] e)^e 2 ^- 3 £(-!)" 



n=l 



l+1 . fnnx \ -s£s±tr, 

n sin e zm^ 

\ M J 



(8) 



The proportionality constant is set by using the normalization constant: 

POO 

dM dt m P(M,t m ,e) = l. 



'x JO 

It is easy to show that the proportionality constant A(e) = e~ 2 . Thus, in the limit 
e — > 0, we finally obtain 



TV 

M 3 



£(-i) n 



71=1 



i /mrx \ n 2 ^ 2 , 

n sin e 2« 2 r 

\ M J 



(9) 
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As a first check, let us show that J °° dt m P(M,t m ) = Xo/M 2 , thus recovering the 
marginal P(M) of the maximum in Eq. HI Integrating over t m , we get 



P(M) = — f sin (^) = *L (10) 

y 1 ttM ^ n V M J M 2 K ' 

n=i 

where the last identity can be found (and derived easily) in [T7] . 

Finally, from Eq. [9] an integration over M (note that M varies from xq to oo) yields 
the desired marginal P(t m ): 

P(t m ) =n — n sin (_) e ^ 



2^- / du cos(u)e 2x o . (11) 



The sum in Eq. [TT] can be expressed in terms of a known special function and we get 

(12) 



P (tm) 



2-Ktr 



where ■d 4 (z,q) is the fourth of Jacobi's Theta functions ([19]). Subsequently one can 
obtain the large and small t m asymptotics of P(t m ) from the exact expression in Eq. [T2j 

Large-t m asymptote: We first consider the case when t m 3> Xq. Changing variables in 

Eq. [TT] through z — x ttiU and letting z — > gives for t m 3> xi: 

Small-t m asymptote: In the opposite limit t m -C x^, we start from Eq. [9] and first take 
a Laplace transform: 



/OO i") °° 

dt m e~ st ™P (M, t m ) = — ^(-l) n+1 
71 n=l 



sm vw) 



sinh(xov 2s) 



[n 

Msinh(Mv / 2s)' 

where the sum of the series can be found in [17J. Letting s become much larger than 
e obtain: 

oo e -VTs(M-x ) 



x Q 2 and M 2 , we obtain: 



poo 

/ dt m e- st ™P(M,t r , 
Jo 



M 

which, after the Laplace inversion ([18]) yields: 

PIM.^g^pV^. (14) 
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Integrating over M gives for t m <C x^: 

1 



(15) 



x y/2nt m 

Thus, P(t m ) has power law behavior at both large and small tails. For large t m , 
the probability density falls off as P(t 

rn) ~ tm ^ 5 whereas for small t m it diverges as 
P{t m ) ~ tm ■ The exact analytical form of P(t m ) and its asymptotes from Eqs. fT2l [T3l 
[15] are plotted (using Mathematica) in Fig. 4 together with the points obtained from the 
numerical simulation (with 1,000,000 realisations). They are in good agreement with 
each other. 
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P(tJ [Eq. 12] 






Small-t m asymptote [Eq. 15] 






Large-t m asymptote [Eq. 13] 
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Figure 4. The probability density P(t m ) and its asymptotes for the driftless case. In 
the simulation, Xq was set to 1. 
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2. In the Presence of a Negative Drift 

We now consider a Brownian motion in the presence of a drift \i. For fi > 0, it is clear 
that the walker will excape to oo with a nonzero probability. This means that with a 
finite probability t m — > oo. Therefore this case is not much of interest in the present 
context. Instead, we focus here on the opposite case where the drift is towards the 
origin, i.e., \x < 0. The Langevin equation describing the motion becomes: 

f = -M + {(*). 

where is the Gaussian white noise with (£(£)) = and (£(£)£(£')) = S(t — t'). We use 
the same strategy as in the driftless case, i.e., splitting the motion into two independant 
parts (before and after t m ) and introducing a small cut-off e. 

On the right-hand side: Letting, as in the driftless case, q(x) be the probability that a 
Brownian particle starting from x G [0, M] exits the interval for the first time through 0, 
we have as before: 



q(x) 



q(x + Ax)4>&t(Ax) dAx. 



(16) 



In the presence of a drift , the mean value of Ax is no longer and one can easily show 
that the analogue of Eq. [2] now reads: 

q"{x) - 2\n\q(x) = 0, q(0) = 1, q(M) = 0. (17) 

The solution is: 

q ( x ) = sinh(|/xKM-x)) ^ |x h ^ 



and so we have: 

q(M 



sinh(|/i|M) 



sinh(|/x|e) 



M(M-s) 



(19) 



' sinh(|/^|M) 

As in the driftless case, the probability density of M can be obtained by differentiation 
of Eq. [18] with respect to M, as was done in [15] : 



V ; sinh 2 (|/i|M) 



(20) 



where we have added the subscript "d" to indicate that the density corresponds to the 
drifted case. 



On the left-hand side: We use the same path integral method as in the driftless case. 
The weight of a path is now proportional to: 



exp 



dx 

(It ( — + 



exp 



\ur i"'"' dx 

-tm - \lA / dr — 



exp 



tm , fdxV 
dT \Tr) 



(21) 
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The position of the Brownian particle at t = and t = t m is known, so we can substitute 



(M — e) — Xo for J Q m dr^p in the first exponential factor on the right-hand side of Eq 
The propagator for the drifted case will therefore be equal to that for the driftless 



case (given in Eq. [5]) multiplied by the factor exp 
be given by: 



fi\x Q -^-t m -\fi\{M-e) 



exp 



\fJL\X 



I At I 



— y^sin 
M ^ 



nnxo 
M 



sm 



nir(M - e) 
M 



and will 



(22) 



As in the driftless case, we multiply the weights of the left and right side of t r 
derived above (Eqs. [T9l and 122]) . and take the e —>■ limit to obtain: 



P d (M,t r 



x - 



-P(M,t m ) 



(23) 



sinh(|/i|M) 

where P(M, t m ) is the joint density for the driftless case given in Eq. (Q. Once again, by 
integrating over t m , one can recover the marginal probability density of the maximum 
P d {M) derived originally in [15]. On the other hand, integrating over M gives the 
marginal Pd(t m ). We were not able to derive a compact expression for Pd(t m ) as in 
the driftless case, though the asymptotes of Pd(t m ) can be derived explicitly as shown 
below. 



Small-t m asymptote: From Eq. [23], we can derive very quickly the behaviour of Pd(t m ) 
when t m <C x^. Substituting in Eq. [23] the asymptotic result for the driftless case from 
Eq. [14] and integrating over M we get: 

P d {tm) ~ \ | | ; /9 - r - (24) 
smh(|/i|x ) y/27Tt m 

Thus for small t m , Pd(t m ) diverges as t m , as in the driftless case. 

Large-t m asymptote: To study the behaviour of Pd{t m ) when t m 3> Xq, we start from 
the following expression for Pd{t m ): 



/"OO 

P d {t m ) = / dMP d (M,t r , 

JX 



l„IW ll»|io- Ji ^« 



XQ smh (fxM) 

2 

00 Iz/IttpH^o- 1 ^-*- x 



rfM IPI - , — r — -> (-l) n+1 nsn -) e ^ tm . (25) 

sinh u M M 2 ^ V ; \ M ) K J 



n=l 



The series in Eq. [25]is dominated by the first term (n = 1) for large t m . Hence, retaining 
only the n = 1 term and making a change of variable y — 1/M in the integral, we get: 

j^o-l^t™ /"^ j.. sin(7rs y) 



/o smh(|/i|/y) 
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For large t m , the most important contribution to the integral comes from the small y 
regime. Expanding the sin and sinh functions and keeping only the leading order term 
reduces the integral to: 



27T 2?o|^|e'' 



dyye 



it ,.2 i M 



ytn 



Letting h(y) = ^y 2 + we next use the saddle point method to obtain the leading 
term via minimizing the function h and get: 



Pd{t r 



,2s .4 

- w — 7T6aT |M| 3 e 



\n\xo 



-5/6 -ML tm -§ (l^f^tU 3 



(26) 



Thus, as expected, the density Pd(t m ) has an exponential decay for large t m in 
presence of a negative drift. Figure 5 shows a plot of the asymptotes (Eq. 1241 and 1261) 
together with the data from numerical simulation (1,000,000 realisations with = 0.1). 
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Figure 5. The simulated probability density P(t m ) and its asymptotes in the presence 
of a drift towards the origin = 0.1). In the simulation, xq was set to 1. 
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3. Summary and Conclusion 

In summary, we have obtained an exact expression for the probability density of the 
time at which a Brownian motion attains its maximum before passing through the 
origin for the first time, and studied the tails of this probability density both for the 
driftless and for the drifted Brownian motion. This was done by first computing the 
joint distribution P(M, t m ) of the maximum M attained and the time t m at which 
it is attained. In the context of the queuing theory, the result that P(t m ) decreases 
monotically with increasing t m suggests that the beginning of a busy period is more 
likely to be the time at which a queue is at its longest. 

It would be interesting to derive the explicit results, obtained here by the path 
integral method, from the general theory of nitrations in Brownian motion developed 
recently in [201 [21]. 

It would also be of interest to extend this calculation to the discrete-time random 
walk case which remains a real challenge. 
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